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ABSTRACT 

RX J0513. 1+0851 and RX J0539. 9+0956 were previously identified as young, 
low-mass, single-lined spectroscopic binary systems and classified as weak-lined 
T Tauri stars at visible wavelengths. Here we present radial velocities, spectral 
types, v sin i values, and flux ratios for the components in these systems resulting 
from two-dimensional cross-correlation analysis. These results are based on high- 
resolution, near-infrared spectroscopy taken with the Keck II telescope to provide 
a first characterization of these systems as double-lined rather than single-lined. 
It applies the power of infrared spectroscopy to the detection of cool secondaries; 
the flux scales as a less steep function of mass in the infrared than in the visible, 
thus enabling an identification of low-mass secondaries. We found that the RX 
J0513. 1+0851 and RX J0539. 9+0956 primary stars are fast rotators, 60 km s -1 
and 80 km s _1 respectively; this introduces extra difficulty in the detection of 
the secondary component as a result of the quite broad absorption lines. To 
date, these are the highest rotational velocities measured for a pre-main sequence 
spectroscopic binary. The orbital parameters and mass ratios were determined by 
combining new visible light spectroscopy with our infrared data for both systems. 
For RX J0513. 1+0851, we derived a period of ~4 days and a mass ratio of q = 
0.46 ± 0.01 and for RX J0539. 9+0956, a period of ~1117 days and a mass ratio 
of q = 0.66 ± 0.01. Based on our derived properties for the stellar components, 
we estimate the luminosities and hence distances to these binaries at 220 pc and 
90 pc. They appear to be significantly closer than previously estimated. 
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Introduction 



Stellar mass is a vital parameter in the life of a star and determines all the evolutionary 
stages of a single object. Theoretical models for the evolution of young stars are widely used 
to estimate stellar masses and ages. For the pre-main sequence (PMS) phase, the theoretical 
calculations for stars above 1.2 M & are approximately consistent with each other. However, 
for stars less massive than 1.2 M G) , the predicted masses are different for different models 



( Hillenbrand fc White! 12004 ; ISimonll2008l ). Thus, for low- mass young stars, these estimates 
are not reliable and few observations exist with which to improve the models. 



Close to half of stars with spectral type GKM have a stellar companion (e.g.. lRaghavan et al 



2010l ). In the last two decades, theoretical formation models for binary stars have been evolv- 
ing. These models assume that the two stars were formed together in the same environment 
and that th eir characte ristics strongly depend on the initial conditions within the molecular 
cloud core (IBatd 120001 ) . The predictions of these models point out that close binaries are 
more likely to have similar mass components. However, it is hard to determine the accuracy 
of these calculations because a direct test requires looking at the statistical properties of 
a large sample of young double stars, such as the mass ratio distribution. Exi sting obser 



vations reveal at least some fraction of binary systems with small mass ratios (IPrato et al. 



2002a]). However, these results are limited by the small number of young binary systems 
with accurately measured mass ratios. 

The study of individual components in young spectroscopic systems offers further insight 
into binary formation. The projected rotational velocity of a star provides clu es to the angu 



lar m omentum evolution and may be correlated with the age of the system (IDelfosse et al. 



19981 ) and the initial conditions of the p arent molecular cloud. Some studies of t he rotational 



veloci ty indicate strong stellar activity (jBrowningj|2008t iMeibom et al.ll2007l ). as lAlcala et al. 
f)l996h suggest with respect to RX J0513.1+0851. 



Young spectroscopic binary (SB) systems are fundamental for the understanding of the 
formation and the evolution of pre-main sequence stars. The dynamical mass ratio can be 
calculated from the radial velocities (RVs) of a double-lined SB system (SB2). For systems 
that are angularly resolvable o r for eclipsing systems, the comp onent stellar masses may be 
measured with high accuracy ( jSchaefer et al.l 12008k l2012t iBoden et al.l |2012| ) . Many of the 
known young SB systems are single-lined (SB1). An SB1 does not provide information about 
the mass ratio; it only allows us to determine the mass function from the calculated orbital 
parameters, based on the observed primary star RVs. The orbital incli nation and masses of 
the primary and secondary stars are inseparable in the mass function ( iMazeh et al.l 120021 ) . 



Historically, most spectroscopic observations of SBs have been made in visible light. In 
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order to increase the sample size of young SB2s with measured mass ratios, we used the fact 
that the low-mass secondary component is cooler and redder than the primary component 



and th us can be detected more efficiently at infrared (IR) wavelengths (e.g. iPrato et al. 
2002bh . Applying this technique, we have observed RX J0513. 1+0851 and RX J0539.9+0956 
spectroscopically in the IR. 



Alcala et all (119961 ) characterized RX J0513.1+0851 and RX J0539.9+0956 as young, 



with ages of < 10 7 years, and with spectral types of K3 and K5, respectively. Selected from a 
survey made by the X— ray satellite ROSAT, RX J0513. 1+0851 showed an emission line with 
equivalent width of W(Ha) = —19.5 ± 0.30 A, an absorption line of W(Li) = 0.25 A, but an 
absence of forbidden lines such as [SII] AA 4068, 4132, 6717 and 6731 A, and [OI] A 6300 and 
6363 A, characteristic of classical T Tauri stars (CTTs). Base d on these res u lts, A lcala et al. 



classified this system as a weak-lined T Tauri star (WTTs). lAlcala et al.l (119961 ) attribute 
the Ha emission as indicative of a very active chromosphere or flare activity during the 
spectroscopic observations. In addition, RX J0513. 1+0851 is not extinguished and does not 
present a near-IR excess (Table [TJ, which indicates a lack of circumstellar or circumbinary 
dust. 



The distance to RX J0539. 9+0956 was estimated at 460 pc (lAlcala et al.lll996[). with 



a tem perature of approximately 4,840 K and a stellar luminosity of 15.9 L Q (IMarilli et al. 



20071 ). Martin (1997) found a very short rotation period (Tabled]) for this system, possi- 
bly symptomatic of an evolved low-mass WTTs approaching the zero age main sequence 
(ZAMS) with a steady increase of the rotation rate during its contraction. T his object was 
treated as a single star to determine its rotational period (jMarilli et al.l 120071 ). a reasonable 
assumption for visible light observations. Table [J summari zes the general properties of RX 
J0539.9+0956 JMarilli et al.ll2007l : lAlcala et al.lll996l l200oh . Near- and mid-IR photometry 
for RX J0539.9+0956, from 2MASS and Wide-field Infrared Survey Explorer (WISE) data, 
indicates a lack of significant IR excess and thus little or no circumstellar or circumbinary 
dust. 

In this paper, we present the detection of low-mass secondary components in two rapidly 
rotating, young spectroscopic binary systems, RX J0513. 1+0851 and RX J0539. 9+0956, 
found by using high-resolution IR spectroscopy to measure the component RVs. These 
systems are unique within the small sample of about 50 known PMS SB2s because of their 
very rapid rotation. Section 2 provides a summary of our observations and data reduction. 
We describe our analysis in section 3 and our results in section 4. Section 5 is a discussion 
and section 6 lists our conclusions. 



-4- 



2. Observations and Data Reduction 

2.1. Visible Light Spectroscopy 

High-resolution visible light spectra for RX J0513. 1+0851 were taken over a period of 6 
years. A total of 32 spectra were obtained between 1996 D ecember and 2002 November with 



two nearly identical echelle spectrographs (ILathamlll992f ) mounted on the 1.5m Tillinghast 



reflector at the F. L. Whipple Observatory on Mt. Hopkins (AZ) and on the 4.5 m-equivalent 
Multiple Mirror Telescope (also on Mt. Hopkins) before its conversion to a monolithic mirror. 
A single echelle order spanning 45 A was recorded with photon-counting Reticon detectors 
at a central wavelength of 5187 A with a resolving power of R = A/AA ~ 35,000. The 
signal-to- noise ratios of these spectra range from 9 to 12, per resolution element of 8.5 km 



s" 1 . 



Visible light spectra for RX J0539. 9+0956 were also taken with the two echelle spec- 
trographs on the 1.5m Tillinghast reflector and on the 4.5m-equivalent Multiple Mirror 
Telescope. Additional spectra were taken with the 1.5m Wyeth reflector at the Oak Ridge 
Observatory (in the town of Harvard, MA). A total of 46 spectra were gathered for this star 
using the three telescopes between 1992 January and 2009 January; signal-to-noise ratios 
ranged from 10 to 29 per resolution element and R = A/AA ~ 35, 000. 



2.2. Near-IR Spectroscopy 

IR spectra of RX J0513. 1+0851 and RX J0539.9+0956 were obt ained with the Keck 



II tel escope and the near-IR spectrometer, NIRSPEC, on Mauna Kea (IMcLean et al.lll998 



20001 ). NIRSPEC is a cross-dispersed cryogenic echelle spectrometer using a 1024 x 1024 
ALADDIN III InSb array detector. The wavelength range in the near-IR covers 0.95 to 5.4 
//m; the resolving power is 30,000 (8 km s _1 ) using a 0.288 arcsec slit width (2 pixels). The 
read noise is ~25 electrons rms. The observations were taken in a pattern of four exposures, 
nodded between two positions on the sky, and hence on the spectrograph slit, in an ABBA 
sequence for each target in order to subtract the background emission and bad pixels. 

For the detection of the secondary component, the IR spectroscopic technique requires 
well-defined absorption lines in the stellar spectra. Blending with the telluric linefl inherent 
in the Earth's atmosphere can confuse the secondary star signal and also requires observations 
of a telluric standard star to remove these features. In order to avoid telluric absorption, the 



Telluric absorption lines produced by the absortion of O3, O2, and H2O are common in IR spectra. 
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targets were observed in the if-band (1.545-1.565 /an) with a central wavelength ~ 1.555 /zm, 
correspondi ng to NIRSPEC orde r 49. This order provides spectra nearly free of terrestrial 
absorption (IRousselot et al.ll2000l ). 



IR Observations for RX J0513. 1+0851 were made during a period of 10 years, between 
2002 February and 2012 January, and for RX J0539.9+0956 between December 2001 and 
January 2012; specific dates appear in Tables H] and [5j All data were reduced using the 
REDSPEC software package designed at UCLA specially for the analysis of NIRSPEC data. 
The reduced and barycentric corrected spectra for ord er 49 are shown in Figures [T] and HJ 
Further details of the data reduction may be found in iMace et al.l ( 120091 ). 



3. Analysis 

3.1. Visible Light RVs 

Primary star RV measurements based on the visible light observations of RX J0513. 1+0851 
and RX J0539.9+095 6 were obtained by cross-correlation using the IRAF@ task XCSAO 
(IKurtz fc Mink! Il998l ) and templates from a large library of synt hetic spectra calculated 



by J. Morse, based on model atmospheres by R. L. Kurucjfl (see Nordstroem et al. 1994 



Latham et al.l |2002| ) . Of the four parameters that characterize these templates (effective 



temperature, rotational velocity, metallicity, surface gravity), the two that affect the RVs 
the most are the temperature and rotational velocity. We determined these by seeking the 
best overall match to the observed spe ctra as measured b y the correlation coefficient aver- 
aged over all exposures of each star (see lTorres et al.ll2002l ). For the unknown metallicity we 
adopted solar composition, and we also held the surface gravity fixed at the value log g = 4.0, 
typical for pre-main sequence stars. For RX J0513. 1+0851, we found that a template with 
an effective temperature of 5000 K and a projected rotational velocity of 50 km s -1 produced 
the highest cross-correlation coefficient. In the case of RX J0539. 9+0956, the best match 
was found for a temperature also of 5000 K, and a projected rotational velocity of 60 km 
s _1 . Tables [2] and [3] show the barycentric Julian dates of the observations and the primary 
RVs and uncertainties found for RX J0513. 1+0851 and RX J0539. 9+0956, respectively. The 
uncertainties were originally determined taking into account both the internal errors as well 
as the signal-to-noise ratio of the individual observations. They were then adjusted through 



2 IRAF is the Image Reduction and Analysis Facility, a general purpose software system for the reduction 
and analysis of astronomical data. 

3 Available at http://cfaku5.harvard.edu. 
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iterations by applying a multiplicative factor until the reduced x 2 from a single-lined orbital 
solution based only on these velocities was near unity. These are the uncertainties listed in 
Tables [2] and EJ 



3.2. Infrared Radial Velocities 



To meas ure the IR RVs (Tables 4 and 5) we used a two-dimensional cross-correlation al- 



gorithm (e.g. JZucker &: Mazehlll994l ). The observed spectrum of an SB2 is correlated against 



a model composed of two template spectra, which represent the primary and the secondary 
components. This requires high quality template spectra with considerable similarity to the 
spectra of the stars in the SB2. For this purpose, observed spec tra of high signal-t o-noise 
ratio were used as primary and secondary star templates (e.g., iPrato et al.l l2002ai ). The 
template spectra were corrected for barycentric motion and the continua flattened. It was 
not possible to determine consistent results on the basis of synthetic template spectra (see 
Mace et al. 2012). 

In addition, to obtain a good match between the template spectra and the observed 
spectrum, we rotatio nally broadened the template spectra by convolving t hem with a line- 
broadening function (Prato 2007 ; Bender et al.l2005 ; Gray! 1992 ; Claret 2000 ). The secondary 
to primary light ratios (a = I2/U ) play an important role in maximizing the correlation 



coefficient (IZucker fc Mazehl Il994l ). Therefore, the light ratio was determined using the 
phase at which the component RVs have the greatest difference, and was then held fixed for 
the analysis of all other phases. 

Uncertainties in the I R RVs for both SB systems were initially determined following 
the same procedure as in iMace et al.l (120121 ); we found primary and secondary star RV 
uncertainties of 4.0 and 5.5 km s _1 and 3.0 and 4.0 km s _1 for RX J0513. 1+0851 and RX 
J0539. 9+0956, respectively. The relatively large uncertainties are likely the result of the high 
rotational broadening in the absorption lines. In the derivation of the double-lined orbital 
solution, we repeated the same procedure as for the visible light data alone and iteratively 
multiplied the uncertainties by multiplicative constants until the \ 2 value was close to unity 
(§3.3). 



3.2.1. RX J 05 13. 1+0851 



The observed template spectra (IPrato et al.ll2002af ) that best matched our RX J0513. 05+0851 
observations were HD 283750 and GL 752 A, of spectral types K2 and M2, respectively. We 
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found projected rotational velocities, v sin is, of 60 km s _1 for the primary component and 
30 km s _1 for the secondary component, and an intensity ratio between the two components 
of 0.20 (Table [6]). These projected rotational velocities and the intensity ratio maximize the 
cross correlation coefficient. We estimated an uncertainty of ± 10 km s _1 for the v sini values 
and ± 0.10 for intensity ratio by varying these parameters until the correlation coefficient 
decreased measurably. Table 0] shows the barycentric Julian dates of the observations, along 
with the RVs and final cr(RV) values fo r RX J0 5 13.05 +0851. As a preliminary analysis, we 
made use of the procedure described in IWilsonl (1194 ll ) to derive approximate values of q = 
0.44 ± 0.01 for the mass ratio and 7 = 68.90 ± 0.52 km s _1 for the center of mass velocity 

(Fig. ED. 



3.2.2. RX J 0539. 9+ 09 56 

From Figure[2]it is evident that the absorption lines of RX J0539. 9+0956 are very broad, 
meaning that there is a high rotational velocity in the dominant stellar component. Cross- 
correlation of the spectra for each epoch yielded a K5 (GL 1094) with v sini = 80 km s _1 
and an M3 (GL 15 A) with v sinz = 40 km s _1 for the primary and secondary, respectively. 
The best fit was found with an M3/K5 intensity ratio of 0.25 (Table [6]). Uncertainties of ± 
10 km s _1 and ± 0.10 for the vsmi values and intensity ratio, respectively, were found as 
for RX J0513. 1+0851. Table [5] shows the IR RVs, final cr(RV) values, and barycentric dates 
of the observa tions for RX J0539. 9+0956. A preliminary analysis was performed following 



Wilson! ( 119411 ); we plotted the six epochs of the primary versus secondary RVs (Figure HJ) 



and found a mass ratio of q = 0.69 ± 0.08 and a center-of-mass velocity 7 =15.7 ± 0.9 km 



s" 1 . 



We note a similar inconsist ency in the magnitu de of the v sin i values between the visible 



light and IR observations as in iMace et al.l ( 120121 ); the IR values are systematically larger 



This difference possibly results from how v sin i is taken into account in the suite of cross- 
correlation templates. For the visible light templates, synthetic spectra were created with 
rotation inherent in the model atmosphere. To simulate a range of vsini values in the 
observed IR templates, a rotation kernel was convolved with the template spectra (e.g., 



Pratol 120071 ). Although it is unclear why this might result in discrepant values for the vsini 



measured in the V- and if -bands, it seems unlikely that the difference has a physical origin. 
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3.3. Results 



We used the Levenberg-Marquardt approach ([Press et al.l Il992l ) to provide a least- 
squares orbital fit on the basis of the measured RVs. We excluded IR spectrum taken 
for RX J0513. 1+0851 on UT 2002 Jan 1 from our analysis because, in spite of the high 
signal-to- noise (Figure 1), the secondary star RV was completely inconsistent with our so- 
lution based on all other epochs. In order to adjust the measured visible light and IR RVs 
so that the reference frame zero points from the two data sets were consistent, we found 
orbital solutions on the basis of the combined data with an extra free parameter, A7, which 
indicates the difference in the center of mass RV between the visible and IR RVs. We then 
corrected all the RVs to a common center of mass velocity by adding this constant, A7, to all 
the measured IR RVs in order to reference the results to the frame of the visible light data. 
For RX J0513. 1+0851 this difference in 7 was -0.5 km s" 1 and for RX J0539.9+0956 it was 
—0.96 km s _1 (Table [7j). Combining the visible light and corrected IR RVs produced the 
best orbital parameters. Table [7] lists the orbital period (P, in days), eccentricity (e), time 
of periastron passage (T ), longitude of periastron (u, in degrees), velocity semi-amplitude 
for the primary and secondary components (K±, K 2 , in km s -1 ), and center of mass velocity 
(7, in km s _1 ). Derived quantities are the projected semi-major axis of the primary and the 
secondary (a 12 sin i), the mass ratio (q = M 2 jMi), and the minimum masses (Mi sin 3 i and 
M 2 sin 3 i). We cannot determine Mi, M 2 or i independently. 

As noted in § 3.1 and § 3.2, the best orbital solutions were obtained with a series of 
iterations, multiplying first the visible light RV uncertainties, then the visible light plus IR 
primary star RV uncertainties, and finally the IR secondary star RV uncertainties by con- 
stants such that the visible-only single-lined solution, the visible+IR single-lined solution, 
and finally the full visible+IR double-lined solution x 2 P er degree of freedom was approx- 
imately unity. Respectively, the constants found for RX J0513. 1+085 were 1.25, 0.98, and 
0.9. For RX J0539. 9+0956 we found 1.2, 0.5, and 0.77. The final, multiplied uncertainties 
used in the non-linear fit are shown in Tables [2j [3j H] and [5j Uncertainties in the orbital 
parameters were determined from the orbital solution found using the Levenberg-Marquardt 
approach. 

The mass ratio and the center of ma ss veloci t y for RX J0513. 1+085 given in Table [7] 



are consistent with those derived from the IWilsonl (1194 ll ) approach to within 1-2(7. The RV 



versus phase curve is plotted in Figure 0, showing the primary star velocities in the visible 
and IR, and secondary star velocities in the IR, as a function of phase. The observed RVs 
are almost all within la of the calculated orbit. 



Simi larly, fo r RX J 0539.9+0956, the mass ratio and the center-of-mass velocity derived 
from the IWilsonl (119411 ) approximation are consistent with those in Table [7J to within 1-2<t. 
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The orbital solution is plotted in Fig. |6] the observed RVs are typically within la of the 
orbital solution except for a few of the visible light RVs. 



Discussion 



4.1. RX J0513.05+0851 

We measured the primary and secondary stars in the RX J05 13. 1+0851 SB by cross- 
correlating a K2 and an M2 with respective projected rotational velocities (v sin i) of ~ 60 km 
s _1 and ~ 30 km s _1 . This factor of two difference in their projected velo cities could result 



from differences in accretion or in magne tic braking as a function of mass (IKonopacky et al. 
20121 ; iGomez Maqueo Chew et al.ll2009l ). It is also possible that the stellar rotation axes 
are not aligned, alt hough this would be surprising in such a short-period (~4 day) system. 
Alcala et al.l ( 119961 ) lo cated this sy stem at ~ 400 pc, corresponding to the distance of the 



Orion cloud complex (Jeffries! 2007 ). Using 2MASS JHK magnitudes to determine an A v of 
~ 0.38, as in lPrato et al.l ( 2003 ). and assuming that RX J0513. 1+0851 is in Orion, the total 
luminosity is L ~ 5.5 L Q . This is very large for these spectral types, even if we suppose the 
system is only a few Myr old. In addition, RX J0513. 05+0851 has a high center of mass 
velocity, 7 ~ 66.8 ± 0.6 km s -1 , which i s not characteristic of Ori on or of the known young, 
nearby stellar associations and clusters ( jZuckerman Songll2004f ). These discrepancies lead 
us the following considerations: 



RXJ0513. 1+0851 is probably not part of the Orion cloud complex. iDolan fc Mathieu 
( 1200 ll ) associated strong lithium sources identified in A Ori with RVs of the systems 
and found a strong peak at ~ 24.5 km s" 1 with a dispersion of 2.3 km s _1 (Figure 
[7]). Previously, the RV ran ge found for some Orion member stars was 21 to 30 km 



s 1 OHartmann et al.lll986l ). With a center of mass velocity of 66.8 ± 0.6 km s 



RX J0513. 1+0851 is many sigma outside this range. RX J0513. 1+0851 is likely closer 
than 400 pc an d part of a differ ent population. We examined the proper motion of RX 
J0513. 1+0851 (IHog et al.ll2000f). and compared it to ^a dozen putative 32 Ori members 
(iMace et al.lbood ; Ishvonski et al.lboid lMamaiekl[2012T ). Although RX J0513.1+0851 
is in the same part of the sky as 32 Ori, its center of mass velocity and proper motion 
are distinct; thus, we conclude that it is not a member of this group either. 



For RX J0513. 1+0851, W(Li) = 0.25 A (lAlcala et al.lll996h . on the lo w end of the range 
of lithium-rich stars in Orion with values between 0.2 and 0.8 A ( Dolan Sz Mathieu 



200 ll ; lAlcala et al.l 120001 ). Figure [8] shows lithium equivalent width versus effective 
temperature for single stars. Using the derived spectral type for the primary star 
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of K2 and assuming the system's W(Li) measurement is associated with the primary 
star, the location of RX J0513. 1+0851 prima ry in Figure [8] is con s istent with the upper 
envelope for young open clusters adopted by lMartm fc Magazzul (119991 ). Traditionally, 
low-mass PMS objects are found in the vicinity of star-forming regions. However, there 
is evidenc e for young Li-ri ch, magnetically active, late-type WTTs which are widely 
dispersed (jNeuhauserlll997l ) . If RX J0513. 1+0851 is a PMS system and part of an open 
cluster, it is likely still moving out from the region where it was formed. Alternatively, 
RX J0513.1+08 51 could have formed locally in small cloud- let, which has dispersed 
jFeigelsonlll9~96L 



Authors such as [Sterzik et al.l ( 119951 ) and iReipurth fc Clarkd ( 120031 ) have suggested a 
scenario of possible ejections of low-mass stars from the molecular parent cloud. This 
is based on the dynamical interactions that occur within small groupings, specifically 
the scenario of three stars which pass close enough that their motions are strongly 
perturbed relative to one and other. The interchange of energy within the three body 
system will cause one of the stars to acquire more kinetic energy, resulting in the 
ejection of the lightest member or least stable member. The total energy is conserved 
and the kinetic energy of the escaping star is released by the formation of a binary 
star system. However, the ejected star with the higher kinetic energy should be single, 
thus this does not satisfactorily explain the high center of mass velocity for the RX 
J0513.1+0851 binary. 



In order to estimate a distance we rely on the models of iBaraffe et al.l (119981 ). shown 
in Figure [9j Although these and other tracks are not totally accurate, and our project 
aims to provide fundamental data to improve these models, we use them here as a 
tool to estimate some of the stellar properties. The reasonableness and consistency of 
our results may also be indicative of the accuracy of the models^. Given the spectral 
types determined for RX J0513. 1+0851 in our cross-corr e lation analysis, K2 and M2, 
we adopt the intermediate temperature scale of iLuhmanl (119991 ) and estimate the cor- 
responding effective temperatures (T eff ). Uncertainties for T eS are ~145 K, equal to 
±1 spectr al subclass. Projectin g the values of log(T eff ) in a vertical line (not plotted) 
across the IBaraffe et al.l (119981 ) tracks in Figure [9] and requiring that the two binary 
components be coeval, we found that the 10 Myr isochrone is most consistent with the 
intensity ratio, 0.2, determined from the two-dimensional cross-correlation, given the 
values of the component T eS . The combined T eff values and the 10 Myr isochrone yield 
a mass of 1.2 ± 0.1 M for the primary star and 0.58 ± 0.1 M & for the secondary star. 



The GAIA mission will return astrometric data which will enable a precise measurement of distance to 
these systems within the next 3—4 years. 
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Fortuitously, these values result in a mass ratio almost identical to our value of 0.46. 
The adopted T e g- along with the 10 Myr isochrone lead to estimated luminosities for the 
components of L\ = 1.12 L & and L 2 = 0.16 L , and a total lum i nosity for the system 
of L tot = 1.28 L & . Inverting the process described in lPrato et al.l (120031 ) for finding the 
luminosity based on the near-IR magnitudes, A„, K-band excess, and distance, we now 
use 2MASS magnitudes, A v = 0.38, and zero K-band excess to estimate a distance to 
RX J0513. 1+0851 of 220pc. 



4.2. RX J0539.9+0956 

Results from our analysis revealed spectral types of K5 and M3 and rapid rotation for 
both the primary component of RX J0539. 9+0956, v sin % ~ 80 km s _1 , and for the secondary, 
v sin i ~ 40 km s -1 . The detection of the secondary component was quite challenging because 
of the high component rotational velocities. The difference between the rotational velocities, 
a factor of two, could be attributed to the initial formation condit ions or to uneven early ac- 
cretion iKonopackv et al.ll2012[ bomez Maqueo Chew et al.ll2009( l. as for RX J0513. 1+0851. 
It is also possible, given the ~ 3 year period of this system, that the stellar inclinations are 
distinct. 



Marilli et al.l ( 120071 ) found a luminosity for RX J0539.9+0956 of L = 15.9 L assuming 



this system is a member of Orion and a mean distance to the Orion complex of 460 pc; 
the relatively large lithium equivalent width (Table [T|) is consistent with the corresponding 
assumption of youth. However, stars of spectral types, K5 and M3, are not massive enough 
to emit such a high luminosity, even if we assume an age of 1 Myr. We thus take into 
account similar considerations and follow a parallel procedure as describe in § 4.1 for RX 
J0513+0851. We draw the following conclusions: 



The center of mass RV of RX J0539. 9+0956 is 15.40 ± 0.33 km s _1 , only moderately 
inconsistent with the R V of the A Ori associati on in Orion of 24.5 km s" 1 with dis- 
persion of 2.3 km s _1 (IDolan fc Mathieul l200ll ) (Figure [7]) and the RV ran g e from 



21 to 30 km s 1 found for some Orion member stars by lHartmann et al.l ( 119861 ). 



Thus we cannot necessarily exclude this system from membership in the Orion re- 
gion given the relatively large scatter in the association members' RVs. A stronger 
argument against membership in Orion is that its ~400 pc distance would necessi- 
tate an unphysically large luminosity, given the spectral types of the stars in RX 
J0539. 9+0956. As in the case of RX J0513. 1+0851, we examined the prop er motion of 



2009; Shvonski et al.l 



2010 



RX J0539.9+0956 (IHog et al.ll2000h compa red to those of 32 Ori members (IMace et al. 



Mamajekll2012l ). In spite of the spatial coincidence with 32 
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Ori, the RX J0539. 9+0956 center of mass velocity and proper motion are distinct thus 
it is unlikely to be a 32 Ori member. 

Even if it is not located in Orion, the presence of a small lithium equivalent width 
in RX J0539.9+0956 indicates youth. Figure j suggests th at RX J0539.9+0956 is 
a post T Tauri star (PTTs), as described by iMartml (119971) . If the transition from 
WTTSs to PTTSs takes place at around 10 Myr of age JMartmlll998h . a PTTs could 
be defined as a low- mass young star that has n ot yet reached t he ZA MS. This points 
to a relatively young age for RX J0539. 9+0956. iBriceno et al.l (119971 ) argued that the 
presence of a detectable lithium equivalent width is still consistent with the zero age 
main sequence (ZAMS). They claimed that the majority of the RASS lithium-rich stars 
found in the general direction of star forming regions are not PM S. Another scenario 
for the presence of lithium was proposed by lAlcala et al.l (119961 ) who posit that for 
late-type stars it is possible to overestimate the W(Li) in low-resolution spectra with 
blending between nearby lines, especially for fast rotators. A third possibility is that 
stellar rotation may affect lithium depletion, h ence the preservation of lithium content 
at high v sin i rates ( iBalachandran et al.l 1201 if ) . 



To estima te the RX J0539.9+ 0956 l uminosity and d istance, we again relied on the 
models of iBaraffe et al.l (119981 ) . From iLuhmanl (119991 ) we found values of T e g of 4365 
K and 3468 K for the spectral types determined by cross-correlation, K5 and M3. The 
isochrone most consistent with these T eff s, coeval binary components, and our observed 
flux ratio, 0.25, has an age of about 30 Myr (Figure [9]). The theoretical models predict 
for the primary star a mass of 0.8 ± 0.1 M & and 0.4 ± 0.1 M for the secondary star. 
The corresponding mass ratio, 0.5, is within ~1<t of our dynamical results. Following 
the same analysis as for RX J0513. 1+0851, we find that the corresponding component 
and thus total luminosities (Li = 0.234 L & , L 2 = 0.07 L , and L to tai = 0.3 L & ) imply 
a distance of 90 pc. 



5. Summary 



We used high-resolution IR spectroscopy in order to measure the Doppler shift in 
the absorption lines of both components in two young SB2s, RX J0513. 1+0851 and RX 
J0539. 9+0956. On th e basis of high-resol ution, visible light spectroscopy these systems were 
initially identified by lAlcala et al.l (119961 ) as SBs and were thought to be located in Orion 
cloud complex (d~400pc). Using two-dimensional cross-correlation, we measured the IR 
RVs for the primary and the secondary of both binaries, combined these values with the 
earlier epochs of visible light data, and found a double-lined orbital solution for each system. 
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The best spectral types found in the cross- correlation of the IR data are a K2 with vsini 
= 60 km s" 1 as the primary and an M2 with vsini = 30 km s _1 as the secondary for RX 
J0513. 1+0851, and a K5 with vsini = 80 km s _1 as the primary and an M3 with vsini 
= 40 km s" 1 as the secondary for RX J0539. 9+0956. To date, these are the most rapidly 
rotating young SBs known. From the orbital solution, we found a mass ratio of q = 0.46 
± 0.02 and period of ~4 days for RX J0513. 1+08 51 and q = 66 ± .11 and a period of 
~1117 days for RX J0539. 9+0956. On basis of the iBaraffe et al.l (119981 ) pre-main sequence 
models and our derived stellar properties, we estimated the luminosities, L = 1.28 L Q and 
L = 0.3 L Q , and distances, ~ 220 and 90 pc, for RX J0513. 1+0851 and RX J0539.9+0956, 
respectively. These estimates are important because to date we do not have an accurate 
method to measure the distance and space motion of these targets. These luminosities are 
significantly lower than would be expected if RX J0513. 1+0851 and RX J0539. 9+0956 were 
located at a distance of >400 pc, but are consistent with the low masses and relatively young 
ages of the component stars. RX J0513. 1+0851 and RX J0539. 9+0956 are significant within 
the sample of known young SB2s because of their high rotational velocities. This provides 
us with fundamental clues about the possible initial conditions in the formation of a binary 
system and the subsequent evolutionary phases. 
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Table 1. General Properties of RX J0513.1+0851 and RX J0539.9+0956 



Property RX J0513+0851 RX J0539+0956 



R.A. a (J2000) 


05:13:05.819 


05:39:56.5 


Decl. a (J2000) 


+08:51:31.44 


+09:56:37 


B h (mag) 


13.80 


11.62 ± 0.1 


V h (mag) 


12.59 ± 0.01 


11.17 ± 0.1 


R b (mag) 


11.88 ± 0.01 




J a (mag) 


10.16 ± 0.02 


9.04 ± 0.03 


H a (mag) 


9.51 ± 0.03 


8.52 ± 0.04 


K & (mag) 


9.35 ± 0.02 


8.37 ± 0.02 


Wl c (mag) 


9.251 ± 0.022 


8.267 ± 0.023 


W2 C (mag) 


9.218 ± 0.020 


8.271 ± 0.020 


W3 C (mag) 


9.177 ± 0.034 


8.273 ± 0.030 


W4 C (mag) 


9.111 ± - 


8.367 ± - 


P rot A (days) 




0.764 


W(Li) e (A) 


0.250 ± 0.001 


0.310 ± 0.001 


W(Ha) e (A) 


-19.5 ± 0.1 


-0.19 ± 0.10 



a 2MASS All-Sky Point Source Catalog. 
b SIMBAD Astronomical Database. 
C WISE All-Sky Data Release. 



^Martini (119971 ). 



Alcala et all fll996[ ). 
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Table 2. Visible Light RVs for RX J0513. 1+0851 and Residual from the Orbital Solution 
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Table 2 — Continued 



BJD 




Vl 


a 


(O-C) 


Orbital Phase 






(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 




2452596.8 


832 


26.45 


7.11 


-1.99 


0.7745 


2452601.8 


819 


76.21 


2.6 


0.08 


0.0185 
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Table 3. Visible Light RVs for RX J0539.9+0956 and Residuals from the Orbital Solution 
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0.6545 


2451487.8112 


23.62 


3.20 


1.72 


0.8427 


2451520.7699 


24.35 


5.66 


3.85 


0.8722 


2451540.7491 


24.47 


2.29 


5.03 


0.8900 
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Table 3 — Continued 



BJD 


Vi 


a (O-C) Orbital Phase 




(km s _1 ) 


(km s _1 ) (km s -1 ) 



2451573 


,5729 


17.25 


2.56 


-0.12 


0.9194 


2451592. 


6500 


16.12 


4.03 


0.11 


0.9365 


2451815. 


.8649 


2.68 


3.10 


-2.80 


0.1362 


2451866. 


.8464 


3.01 


2.11 


-3.20 


0.1819 


2452699. 


.5888 


25.60 


2.88 


8.83 


0.9271 


2452937. 


8303 


4.03 


3.38 


-1.48 


0.1403 


2452958. 


.8358 


5.10 


3.11 


-0.64 


0.1591 


2453354. 


.7803 


21.52 


3.78 


2.50 


0.5135 


2453402. 


6936 


20.16 


3.38 


-0.26 


0.5564 


2453425. 


.6159 


17.63 


2.29 


-3.41 


0.5769 


2453448. 


.5899 


19.31 


2.80 


-2.29 


0.5974 


2454043. 


.9197 


7.42 


1.38 


1.96 


0.1302 


2454160. 


.7274 


7.05 


1.76 


-0.82 


0.2348 


2454193. 


6679 


9.85 


1.87 


0.84 


0.2642 


2454514. 


.8194 


20.18 


3.54 


-0.10 


0.5517 


2454845. 


.8486 


23.83 


2.54 


2.14 


0.8479 



-23 - 



Table 4. IR RVs for RX J0513. 1+0851 



BJD 


vi a 


(O-C) 


v 2 b 


(O-C) 


Orbital Phase 




(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 




2452311.75855 


35.8 


4.00 


147.5 


3.04 


0.8177 


2452629.91432 


70.5 


0.43 


65.7 


4.29 


0.9947 


2452678.86556 


105.7 


0.25 


-16.8 


-1.43 


0.1768 


2453032.74867 


108.9 


2.99 


-15.3 


0.84 


0.2450 


2453365.95440 


106.9 


2.10 


-9.4 


4.54 


0.1673 


2455937.93153 


98.3 


-7.90 


-15.9 


1.09 


0.2354 


2455941.89231 


98.6 


-7.93 


-20.1 


-2.40 


0.2211 


a Primary star RV uncertainty is 3.9 km s 1 . 






b Secondary star RV uncertainty is 5.0 km s _1 . 








Table 5. 


IR RVs for RX J0539.9+0956 




BJD 


Vl a 


(O-C) 




(O-C) 


Orbital Phase 




(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 




2452274.92156 


21.1 


0.01 


7.6 


-1.62 


0.5471 


2452311.77812 


22.6 


0.51 


5.9 


-1.82 


0.5801 


2452622.99356 


22.2 


0.04 


6.2 


-1.41 


0.8586 


2453030.94592 


10.1 


1.67 


23.8 


-4.50 


0.2237 


2453365.99928 


21.6 


1.28 


7.3 


-3.08 


0.5235 


2455937.92445 


22.8 


-0.69 


6.0 


0.39 


0.8253 



a Primary star RV uncertainty is 1.5 km s 1 
b Secondary star RV uncertainty is 3.1 km s 
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Table 6. Cross-Correlation Derived Properties 



Object 


Spec. Type 


r~p a 
J-efi 


v sini 


Flux Ratio 






(K) 


(km s _1 ) 




RX J0513.1+0851 


K2 b 


5010 ± 145 


60 ± 10 


0.20 ± 0.10 




M2 C 


3630 ± 145 


30 ± 5 




RX J0539.9+0956 


K5 b 


4365 ± 145 


80 ± 10 


0.25 ± 0.10 




M3 C 


3468 ± 145 


40 ± 5 





a Spectral Type - T e g conversion based on iLuhmanl (119991 ) . 
b Primary Component. 
c Secondary Component. 
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Table 7. Spectroscopic Orbital Solutions from Combined Visible Light and IR Data for 

RX J0513.1+0851 and RX J0539.9+0956 



Orbital parameter 


RX J0513. 1+0851 


RX J0539.9+0956 


P (days) 


4.01829 ± 0.00006 


1117.39 ± 5.34 


T (BJD) 


2452525.46 ± 0.16 


2451663.6 ± 57.0 


e 


0.06 ± 0.02 


0.20 ± 0.05 


K 1 (km s- 1 ) 


39.0 ± 0.7 


9.2 ± 0.7 


K 2 (km s- 1 ) 


84.6 ± 2.7 


13.8 ± 2.1 


u (deg) 


275.8 ± 14.5 


115 ± 18 


7 (km s _1 ) 


66.8 ± 0.6 


15.40 ± 0.33 


A7 a (km s _1 ) 


-0.5 ± 1.5 


-0.96 ± 0.71 


ai sin « (10 6 km) 


2.15 ± 0.04 


138.0 ± 9.6 


a 2 sin i (10 6 km) 


4.66 ± 0.15 


208.0 ± 31.6 


q = M 2 /M 1 


0.46 ± 0.02 


0.66 ± 0.11 


Mi sin 3 i (M Q ) 


0.53 ± 0.04 


0.79 ± 0.27 


M 2 sin 3 i (M ) 


0.25 ± 0.01 


0.53 ± 0.12 



a Factor to adjust visible and IR RVs to the same reference 
frame (§3.3). 
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Fig. 1. — Eight epochs of spectra for RX J0513. 1+0851 in order 49, with barycentric cor- 
rections applied; UT dates of the observations are indicated. An arbitrary additive constant 
was used to offset the spectra for display. 
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Fig. 2. — Six epochs of spectra for RX J0539. 9+0956 in order 49, with barycentric corrections 
applied; UT dates of the observations are indicated. An arbitrary additive constant was used 
to offset the spectra for display. 
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Fig. 3. — Linear fit to t he primary vs. secondary RVs for RX J0513. 1+0851. Results given 
for 7 and q are from the IWilsonl ( ]194ll ) analysis. 
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Fig. 4. — Linear fit of th e prima r y vs. secondary RVs for RX J0539. 9+0956. Results shown 
for 7 and q are from the IWilsonl ( ]194ll ) analysis. 
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Fig. 5. — Radial velocity as a function of phase for RX J0513. 1+0851. The black open circles 
represent the primary star visible light data, the green circles the primary star IR data, and 
the red circles represent the secondary star IR data. The best orbital fit is shown as a solid 
line for the primary and secondary components. The dashed horizontal line represents the 
system's center-of-mass velocity. Uncertainties in the primary IR RVs are 3.9 km s _1 and in 
the secondary star RVs are 5 km s" 1 . 



31 



30 - 



E 



U 
_o 

CD 
> 



O 

a: 



20 - 



10 




J_ 



_L 



0.0 0.2 0.4 0.6 

Orbital Phase 



OA 



1.0 



Fig. 6. — Radial velocity as a funcion of phase for RX J0539.9 +0956. The black open circles 
represent the primary star visible data, the green circles the primary star IR data and the 
red circles represent the secondary star IR data. The best orbital fit is shown as a solid 
line for the primary and secondary components. The dashed horizontal line represents the 
system's center-of-mass velocity. Uncertainties in the primary RVs are 1.5 km s _1 and the 
secondary star RVs are 3.1 km s -1 . 
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Fig. 7. — Figure 4 from iDolan Mathieul ( 1200 ll ) showing the discrimination between field 
stars and low-mass members of the A Ori association. The dashed line indicates an empirical 
cutoff of W\{Li) > 0.2 A, identifyi n g you ng stars with strong lithium absorption. The small 
circles are from IDolan fc Mathieul (120011 ); open symbols are binaries. The green circle and 
red square represent the center of mass RVs for RX J0513. 05+0851 and RX J0539. 9+0956, 
respectively. These symbol sizes are at least 10 times larger than the uncertainties of these 
points (see Table [7j). 
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Fig. 8. — Lithium equivalent width versus effective temperature for sin gle stars. The thick 
and d otted lines represent the upper envelope for young open clusters (IMartm fc Magazzu 
1999). The dashed line indicates the WTTs and PTTs regions as described in Martin 



(jl997 ). The green circle and red square represent RX J0513.05 +0851 and RX J0539. 9+0956 
respectively. 
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Fig. 9. — A Hertzsprung-Russell diagram with the locations of the RX J0513. 05+0851 and 
RX J0539. 9 +0956 comp o nents (§4) compared to the solar metallicity theoretical evolutionary 
models of iBaraffe et al.l (119981 ) . T eff s and luminosities for each component are discussed in 
the text. The evolutionary tracks are shown as dashed lines with mass in units of M . For 
M ^ O.6M the adopted mixing length is a m i X = 1 corresponding to a helium mass fraction 
of Y= 0.275, and for M > 0.6 M the mixing length is a mix = 1.9 and Y = 0.282. Solid lines 
are isochrones representing ages of 1, 2, 3, 5, 10, 30, and 100 Myr (from top to bottom). The 
green circles and red squares represent primary and secondary stars of RX J0513. 05+0851 
and RX J0539.9+0956, respectively. 



